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We present a screened exact-exchange (SXX) method for the efficient and accurate calculation of 
the optical properties of solids, where the screening is achieved through the zero-wavevector limit of 
the inverse dielectric function. The SXX approach can be viewed as a simplification of the Bethe- 
Salpeter equation (BSE) or, in the context of time-dependent density-functional theory, as a first 
step towards a new class of hybrid functionals for the optical properties of solids. SXX performs 
well for bound excitons and continuum spectra in both small-gap semiconductors and large-gap 
insulators, with a computational cost much lower than that of the BSE. 

PACS numbers: 31.15.ee, 71.15.Qe, 71.35.Cc, 78.20.Bh 


The Bethe-Salpeter equation (BSE) [l|,[2| is considered 
the gold standard for calculating the optical properties of 
periodic solids and many other materials. The nonempir- 
ical nature of the BSE guarantees its wide applicability 
and high degree of accuracy, but its computational cost 
becomes prohibitive beyond the simplest systems. Time- 
dependent density-functional theory (TDDFT) M is 
computationally much less expensive, and is therefore a 
popular alternative for the calculation of optical proper¬ 
ties. TDDFT calculations can be orders of magnitude 
faster than the BSE, but none of the existing empiri¬ 
cal or nonempirical exchange-correlation (xc) kernels for 
solids can achieve the same level of accuracy for 

both small-gap and wide-gap solids. The exception is the 
so-called “nanoquanta” xc kernel (2, 11-14], which is as 
accurate as the BSE, but equally expensive. 

Recent TDDFT studies for solids have identified the 
crucial importance of the long-range part of the xc ker¬ 
nel 15,13]. Exact-exchange TDDFT 17. Ijj successfully 
produces excitonic properties, but the Coulomb singu¬ 
larity needs to be cut off, which is equivalent to using a 
screened Coulomb interaction [Till . Hybrid xc functionals 
are defined as a mixture of semilocal (gradient-corrected) 
xc functionals with a fraction of nonlocal exact exchange. 
The B3LYP hybrid functional [2Cj] has been used to cal¬ 
culate optical spectra in semiconductors [U JUl, with 
a generally good description of optical gaps, despite the 
fact that the 0.2 mixing parameter of B3LYP is optimized 
for finite systems. The HSE functional 23jl 24 [ uses exact 
exchange for short-range interactions only; this produces 
very good quasiparticle gaps 25|-27|, but HSE cannot 
yield bound excitons, although it may still give decent 
continuum spectra 28]. 


In this paper we propose a simple, nonempirical 
and material-dependent way of screening the long-range 
Coulomb exchange, which can be viewed as a simpli¬ 
fied BSE approach. We show that this screened exact- 
exchange (SXX) approach outperforms all TDDFT ap¬ 


proaches currently on the market, retaining most of the 
accuracy of the BSE over a wide range of materials, but 
at a much lower computational cost. This builds a bridge 
between TDDFT and many-body theories, and opens up 
new directions towards the development of hybrid func¬ 
tionals for the optical properties of insulating solids. 

Although TDDFT and BSE are very different theories, 
the excitation spectra in solids are in both cases obtained 
through an eigenvalue equation [29j (atomic units [e = 
ft = m e = l/47reo = 1] are used throughout): 
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where i and m denote occupied bands, j and n denote un¬ 
occupied bands, the e’s are single-particle energies (either 
quasiparticle or Kohn-Sham), and ui is the excitation fre¬ 
quency. The main difference lies in the coupling matrix 
Ehxc = Eh + F xc . For optical properties, only vertical 
excitations are considered, so that k, = kj and k m = k„ 
in Eq. ©■ The Hartree part of the coupling matrix is 
the same in the two methods, and is given by 
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The long-range part (G = 0) of the Coulomb interaction 
is omitted so that the eigenvalues of Eq. m correspond 
to poles in the macroscopic dielectric function j2|, l30(. 

For the BSE, as well as for our SXX method, the xc 
part of the coupling matrix can be written as 
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Here, g G G '(q) = -47T7<5 GG ,/|q + G| 2 for SXX (7 is a 
screening parameter, to be further specified below), and 
SGG'(q) = — 4 /TEQp, (q, w = 0) / 1q + G'| 2 for the BSE. 
e ^ 1 is the inverse dielectric function, obtained within the 
random phase approximation (RPA) as 


C GG 


,(q, uj) = Sgg 1 + 


47T 


G| 2 


X ^(q, W ), (4) 


with the RPA response function defined as y RPA = yn + 
XoiA RPA (Xo is the quasiparticle response function |5(). 
In TDDFT, the xc part of the coupling matrix is 
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where / xc ( r,r') = 5v xc (r)/Sn(r') is the adiabatic xc ker¬ 
nel of TDDFT. The structure of Eq. (O is similar to Eq. 
ED, but different from Eq. ED: only the q = 0 part of 
/ xc enters in the expression, so the head (G = G' = 0) 
at q = 0 of the xc kernel plays a much more important 
role in TDDFT than goo(0) in the BSE. 

To illustrate the difficulty in developing a universally 
applicable nonempirical xc kernel, we consider the long- 
range corrected (LRC) kernel @, O, El] 

/xc R GG'(q) = — |q_|_Q|2 ^ GG '’ (6) 


which represents the long-range part of the exact xc ker¬ 
nel in insulators. The empirical parameter a acts as a 
rough approximation to the dielectric screening effects 
within the BSE, but has no clear justification in the 
TDDFT framework, besides giving the correct asymp¬ 
totic behavior. In Ref. [7], the relation a = 4.615 e ” 1 — 
0.213 was proposed, which works quite well for the opti¬ 
cal spectra of semiconductors, but fails for insulators. 

In Ref. f32j, a was fitted against experimental exciton 
binding energies for various materials. It was found that 
the value of a spans a wide numerical range, from 0.595 
(GaAs) to 96.5 (solid Ne). For small-gap materials, the 
relative change in the exciton binding energy caused by 
a change in a is substantial; for large-gap materials, the 
exciton binding energies are not as sensitive. This shows 
how difficult it is to develop a widely applicable nonem¬ 
pirical xc kernel within the TDDFT framework [foil . 

The situation is different in SXX. Unscreened time- 
dependent Hartree-Fock (TDHF) always overbinds exci- 
tons, so 7 has to be in the [ 0 , 1 ] range for the correction 
to be in the right direction. Therefore, it is a much easier 
task to develop nonempirical approximations for 7 than 
for the TDDFT parameter a. 

To derive the SXX screening parameter 7 , we start 
from the self-energy X: 
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where G is the quasiparticle Green’s function, e 1 in Eq. 
CD is the full dielectric screening of the BSE, and we 
want to find a way to average its effect and motivate re¬ 
placing it with a constant. Assuming that we can replace 
the static nonlocal screening e _1 (r, r') with a frequency- 
dependent uniform screening e“ ni (r — r/ , w ), we write 

E(r,r',w) = ^-G(r,r',w) J d 3 r" e“ n 1 i (r-r",u;)u(r"-r'), 

(8) 

which defines the function G. Combining Eqs. CD and 
© in reciprocal space leads to 
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Eq. J9| holds for any q, G, G'. Setting these to zero and 
assuming the functions G and G to be real, we have 
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In the q —> 0 limit, Eq. USD becomes 
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For G/0, Eq. jTTJ) behaves like 0(l/q), so it is impos¬ 
sible to derive a general uniform screening as an average 
of the nonlocal screening of the BSE. However, the head 
(G = G' = 0) of Eq. (fill) has the correct q —>■ 0 behav¬ 
ior, and can be used to approximate the BSE. Assuming 
that G = G, we obtain 


e uni, 0 (q, w ) = e oo 1 (q,°)- (I 2 ) 

Now we make a rather drastic approximation by set¬ 
ting q = 0 , since the long-range interaction dominates 
in solids. The screening parameter 7 then becomes 

7 = e o “o 1 (0,0), (13) 


which is also the inverse of the infinite-frequency dielec¬ 
tric constant, e^ 1 , since phonon effects are not included. 
A similar simplified screening was proposed in Ref. [Uj 
for the nonlocal exchange part of a hybrid xc functional 
in order to obtain good band structures; by contrast, 
the purpose of our SXX is to yield good optical proper¬ 
ties. We use the RPA for eg O L (0,0) in actual calculations. 
Since e = 1 — vx, and the static \ at zero wave vector is 
negative, 7 of Eq. (H3l) is bounded in the [0,1] range. 

Let us now compare the SXX approach with the BSE. 
The main difference is that in BSE the exchange is 
screened by the full inverse dielectric function e^ 1 , which 
makes it much more costly than SXX, where the screen¬ 
ing parameter 7 is just a constant. In practice, a BSE 
calculation is a four-step procedure: (i) ground-state cal¬ 
culation with a diagonalization over the selected fc-points 
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Experimental E b (eV) 

FIG. 1. (Color online) Comparison of calculated and experi¬ 
mental exciton binding energies Eb in various materials (see 
Table Ufor further details). The solid line indicates where the 
calculated and experimental values of Eb coincide. 



FIG. 2. (Color online) SXX screening parameter 7 [Eq. (1131 ] 
and experimental value of e ^ 1 versus the fitted 7 reproducing 
the first exciton, for various materials. B3LYP corresponds 
to a constant value of 7 = 0 . 2 . 


grid (often shifted [35] for optical properties); (ii) quasi¬ 
particle correction, typically within the GW approxima¬ 
tion jUj, but substituted with a simple scissor correction 
step in this work since we focus on the excitonic effects; 
(iii) generation of static screening eQQ,(q,u = 0 ) within 
the RPA; (iv) construction and diagonalization of the 
excitonic Hamiltonian [i.e., Eq. JT|)] containing the in¬ 
gredients listed above. 

Regarding the computational workload, even though 
step (iv) has the worst scaling, step (iii) is often the most 
cumbersome part of the whole procedure, especially when 
one is interested in the small energy region of the spec¬ 
trum and uses the scissor operator procedure: the num¬ 
ber of q-vectors in the screening is proportional to the 
number of fc-points (since q = k — k') even for optical 
properties. This can become very demanding when many 
^-points have to be used together with many G-vectors, 
as is the case for lower-dimensional systems, in partic¬ 
ular 2D [ 37 I ]. In addition, the numerical evaluation of 
e GG' (d- 0 ) has a much worse convergence with the empty 
bands than the evaluation of the spectrum (for instance, 
the screening for LiF requires 20 empty bands, while the 
first exciton peak requires only one empty band). Our 
SXX approach bypasses this step and thus avoids a se¬ 
vere computational bottleneck in the description of opti¬ 
cal properties at BSE-level for complex materials. 

We have calculated the exciton binding energies Eb 
of various semiconductors and insulators with SXX and 
other methods; the results are collected in Table U and 
in Fig. [I] Clearly, SXX produces a much better overall 
agreement with experiment than all TDDFT methods 
(except for GaAs), and yields an accuracy that is com¬ 
parable to BSE across the board. 

The calculations of Eb were done with the method de¬ 


scribed in Ref. [Hj]. We use a scissor-corrected LDA 
ground state (calculated with ABINIT Q) as starting 
point. All calculations use three valence bands and one 
conduction band for Eq. 0. which is sufficient when 
only the exciton binding energies are of interest (but is, 
generally, insufficient for the continuum part). We use 
an 18 x 18 x 18 Monkhorst-Pack grid [39j] for GaAs and 
/3-GaN, a 15 x 15 x 15 grid for MgO, a 10 x 10 x 10 
grid for Ar, Ne, and LiF, and a 20 x 20 x 20 grid for 
other materials. To speed up the calculation, we only 
use the head of the xc kernel when calculating the cou¬ 
pling matrices [i.e., we neglect local-field effects by not 
taking the G and G' sums in Eq. m and ©]. Including 
local-field effects generally changes Eb very little (at most 
-10%). Since the Eb are already small numbers, the re¬ 
sults without local-field effects are sufficiently accurate 
for our purposes. For the calculation of £^(0, 0), we use 
60 bands for GaAs, /3-GaN and MgO, and 30 bands for 
all other materials. 59 G-vectors are used for e -1 , and 
the error for the convergence of e^lO, 0 ) is less than 1 %. 

Figure [2] compares 7 from Eq. m with values of 7 
fitted to reproduce the lowest experimental exciton bind¬ 
ing energies. Aside from a few outliers (such as GaAs), 
the calculated and the fitted values of 7 are very close, 
which explains the good performance for exciton bind¬ 
ing energies in Table [U The experimental e ^ 1 values are 
also plotted in Fig. 0 showing that 6 q O 1 (0, 0) at the RPA 
level is already a good approximation to e z}. It should 
be noticed that the B3LYP hybrid kernel [20j (only the 
long-range part, which corresponds to 7 = 0 . 2 , since the 
calculation only uses the head of the xc kernel) also per¬ 
forms well for semiconductors: Fig. 0shows that 7 = 0.2 
is roughly the average of the semiconductor screening 
parameters. The B3LYP functional was designed with 























4 



GaAs 

/3-GaN 

a-GaN 

CdS 

CdSe 

Ar 

Ne 

LiF 

AIN 

ZnO 

MgO 

Exp. 

3.27 

26.0 

20.4 

28.0 

15.0 

1.90 x 10 3 

4.08 x 10 3 

1.6 x 10 3 

75 

60 

80 

BSE 

— 

— 

172 

66.0 

16.2 

2.07 x 10 3 

3.32 x 10 3 

2.51 x 10 3 

552 

208 

546 

TDHF 

497 : 

1.99 x 10 3 

2.00 x 10 3 

1.28 x 10 : 

* 879 

3.27 x 10 3 

4.26 x 10 3 

4.68 x 10 3 2.37 x 10 3 

1.84 x 10 3 

3.04 x 10 3 

B3LYP2 

0.792 

7.71 

57.4 

48.1 

25.8 

100 

197 

180 

89.4 

55.8 

31.9 

SXX 

0.151 

4.08 

16.4 

26.0 

11.9 

1.33 x 10 3 

3.08 x 10 3 

1.46 x 10 3 

39.9 

30.8 

165 

LRC^ 

0.858 

0.514 

0 

0.513 

1.40 

0.304 

0.127 

1.14 

0 

0.810 

0.076 

B 00 C 

0.332 

0.199 

0 

0.461 

0.895 

1.70^ 

852^ 

32.2^ 

0 

1.09 

0.051 

JG1VP 

0.833 

0.382 

0 

0.741 

1.42 

41.0 

0.593 

993 

0 

4.45 

1.79 


TABLE I. Exciton binding energies calculated with Eq. ©. compared with experimental values (all numbers in meV). All 
calculations are head-only; see text for other technical details. The BSE results for GaAs and /3-GaN were not calculated. The 
estimated error due to the head-only approximation is < 10% for all many-body calculations, and < 5% for TDDFT. 

a Head-only calculation, equivalent to SXX with 7 = 0.2 independent of the material. 
b With the empirical formula of Ref. [?;]. 
c The bootstrap kernel of Ref. [ 8 j. 

d The convergence of the bootstrap kernel strongly depends on the number of bands included in the iterative calculation of the kernel. 

These results are obtained by evaluating the bootstrap kernel with 30 bands. The results reported in Ref. [iidl were not fully converged. 
e The jellium-with-gap model of Ref. T)]. 



FIG. 3. (Color online) Absorption spectrum of LiF calculated 
with SXX and RPA, compared with experiment [40], 



FIG. 4. (Color online) Absorption spectrum of AIN calculated 
with SXX, RPA, and BSE, compared with experiment [Tlj. 
The BSE spectrum of Benedict et al. [|2] is also shown. 


small molecules in mind, so its good performance for 
bound excitons in semiconductors seems fortuitous. 

To demonstrate that our method yields good results 
not only for the exciton binding energies, we present the 
optical spectra of LiF (Fig. [3]), AIN (Fig. |4}, and Si (Fig. 
ED- The spectra are obtained in a standard manner via 
the imaginary part of the macroscopic dielectric function 
(2]. We use 20 bands and 256 fc-points for LiF, 10 bands 
and 256 fc-points for AIN and Si. All calculations include 
local field effects. We obtain a very good agreement of 
the position and strength of the strong bound-exciton 
peak in LiF compared to experiment, which is also ev¬ 
ident from the good agreement between the calculated 


and fitted screening parameters shown in Fig. [2] For the 
smaller-gap materials AIN and Si, the excitonic enhance¬ 
ment of the band-edge spectrum is somewhat underesti¬ 
mated (the bound excitons are not shown in Figs. [4]and[5] 
since Ef, is smaller than the frequency resolution). Here, 
the excitonic enhancement effects in the BSE and experi¬ 
mental spectra are due to continuum excitons. Compared 
to RPA, the SXX spectra in Figs. 0] and [5] give a much 
better description of the excitonic enhancement effects. 

In conclusion, we propose a very simple nonempirical 
screening factor for nonlocal exchange, derived as a sim¬ 
plification of BSE. We show that it is easier to derive a 
good approximation in the many-body framework than 
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FIG. 5. (Color online) Absorption spectrum of Si calculated 
with SXX, RPA, and BSE, compared with experiment [43| . 

developing a better long-ranged xc kernel for TDDFT. 
Our SXX approach yields exciton binding energies of 
a wide range of semiconductors and insulators in good 
agreement with experiment; the performance is consis¬ 
tently better than currently available TDDFT methods. 
The SXX method works well for the optical spectra of 
wide-gap materials, and captures continuum excitonic 
effects in small-gap materials to some extent, although 
there is still some room for improvement. 

The SXX approach constitutes a first step towards a 
hybrid xc kernel specifically designed for optical proper¬ 
ties in periodic insulators and semiconductors. In this 
paper we have focused on the long-range behavior of the 
xc kernel; the next step will be to match the SXX ap¬ 
proach with suitable xc functionals for the short range to 
capture local-field effects. This should have minor effects 
on strongly bound excitons, but is likely to lead to an im¬ 
provement of the continuum part of the optical spectrum. 
Work along these lines is in progress. 
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